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1 Introduction 



The interest of the configuration space structure of a hadronic source has appeared 
primarily in intensity interferometry [Q]. The technique was developed originally to 
estimate the dimension of distant astronomical objects. Since that time it has seen 
widespread application in subatomic physics, in particular in analysis of elementary 
particle collisions f^. The standard HBT procedure involves introducing an ansatz 
describing the geometry of the particle source, on the basis of a physical model. Then 
it investigates its multiparticle characteristcs in momentum space. Many models of 
hadron production based on some phenomenological or theoretical constraints 
have been considered; however, the question as to what is the configuration space 
structure of the QCD cascade when derived explicitly from the fundamental theory 
has not been addressed so far. 

Recently, several groups have analysed in great detail the multiparton distri- 
butions in the QCD gluonic cascades 0. The results of their investigations show 
that perturbative QCD provides a powerful framework not only for the descrip- 
tion of hard quark and gluon jets but also of much softer multiparticle phenomena. 
Although not understood theoretically, the hypothesis of parton-hadron duality p 
provides an apparently successful link between theoretical parton distributions and 
observed particle spectra. This prescription was extensively tested in single particle 
spectra ( and total multiplicities) and found to be in a good agreement with the 
available data (see e. g. 0). Recently, there have appeared indications that it may 
also work for multiparticle correlations These unquestionable successes invite 
one to study further consequences of the theory for processes of particle produc- 
tion. At this point we would like to notice that, if one wants to exploit fully the 
quantum-mechanical aspects of the QCD cascade, it is necessary to study at first 
the multiparton density matrix. The multiparticle distributions calculated so far 
give only diagonal terms of the density matrix and thus represent a rather restricted 
( although very important ) part of the information available from the theory. 

It is perhaps important to stress that in contrast to what is usually believed, 
the interest in studying the multiparticle density matrix is not purely academic. As 
suggested in the density matrix allows one to obtain the multiparticle Wigner 
functions and consequently gives information about the space-time structure of the 
system. This allows one to make predictions about the shape and range of the HBT 
interference with clear experimental consequences. 

In my previous paper I investigated the multiparticle density matrix (DM) of 
the gluonic cascade produced in e~^e~ collision in the framework of double logarith- 
mic (DLA) approximation |jlOl of QCD. I proposed a generating functional to 
obtain integral equations for the multiparticle density matrix in the quasi- diagonal 
limit, i. e. if the energies and emission angles of particles are close to each other. 
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Here I will be presenting the technique for extracting physical information from 
the density matrix approach, deriving the explicit form of inclusive single particle 
density in configuration space. For the sake of simplicity, I will be restricting myself 
to a discussion of its dependence on the transverse space coordinate xt- In solving 
the problem, first I will discuss the derivation of the inclusive single particle density 
matrix d^p{k',k). I will prove that the terms in d^p{k',k) which do not vanish for 
k = k' are leading in the double logaritmic perturbative expansion when as ^ 0, 
momentum of quark ( antiquark ) P generating the gluonic cascade is large : 
P ^ oo and the double logarithmic contributions of the form In P = const 
dominate. Taking into account these leading terms will allow me to obtain the 
explicit analytic form of the density matrix, and consequently the explicit form of 
the single particle density in configuration space. Finally, I will be analysing the 
physical properties of inclusive single particle density, and I hope to find them in 
agreement with intuitive physical expectations. 

2 Density matrix formalism 
2.1 Definition of the density matrix 

In this paper I will concentrate on the single particle density in configuration space. 
I will be discussing it in relation to the QCD-parton cascade within the framework 
of double logarithmic approximation (DLA), using the density matrix formalism 
presented in ||^. Hence in this section I will recall the definition of the density 
matrix for a particle production process, and show the relation between the density 
matrix and particle density in configuration space. Our assumption is that all the 
produced particles are real, i.e. they are on a mass shell. If the production of 
m particles can be realized in different ways represented by a sample of Feynman 
diagrams, then the exclusive m-particle density matrix equals the product of total 
production amplitude S{ki, . ., k^) and its complex conjugate S*(^k'^, . ., fc^) as : 

d (fc^, . ., k^; fci, . ., km) = S {k^, . ., k^)S{ki, . ., km)', (1) 

where the total amplitude S{ki, . ., k^,) is the sum of all contributions S(^£)) from 
Feynman diagrams (D) multiplied by phase space factors (2a;fc)~^/^ : 

m 

S{k^, ..,km) = Yl SiD){h, ..,km) n(2^feJ"'/'- (2) 

D i=l 

For inclusive analysis one constructs the m-particle density matrix as a series of 
integrated n-particle exclusive densities in the form : 

d^'^lk' k' ■ ki k — 
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7 ~ rr [dk]m+i...n c^^^lk'j^, • k^, k^+i, . ., k„; ki, . ., k„); (3) 

n=m ^J- 

where [(iA;]i,,.j = d^ki. .d^kj. This construction scheme imphes, that the diagonal 
elements of the density matrix are equal to particle densities in momentum space. 
There is also an obvious relation between the density matrix and particle density 
in real space. Remembering, that the space-time multiparticle amplitude is the 
on-mass-shell Fourier transform of the momentum amplitude, i.e. 

S{x,,..,xJ = y'[rffc]i...^e*'^'=i*i-^'^^^\.e*'^'='"*'"-*'^-^-5(A;i,..,U; (4) 

where denotes the energy of the ith produced particle, for the multiparticle 
density . .,Xm) we get the relation : 

P ^ (Xi, . ., Xm) = /r,\3rn ^ [c^^] l...m [c^^ ] l...m ^ {ki, ■ ■ , kj^] ki, . . , km) (5) 



(2vr)= 



where the factor ^^^^^3^ was introduced to get the proper normahzation : 

Mx]i...^p'=^/^"(xi,..,a:^) = /[rfA;]i.„^/^/^"(fci,..,fc„). (6) 



2.2 Parton cascade in the double logarithmic approximation 

The double logarithmic approach in momentum space IQ,!!^] gives a good quali- 
tative description of the structure of the gluonic cascade. It accounts only for the 
leading DL contributions to the multiparticle cross section. Although emitted soft 
gluons violate the energy and momentum conservation rules, however, at high ener- 
gies the approximation reproduces quite well the selfsimilar structure of the gluon 
radiation. Let us consider in the framework of DLA the gluonic cascade generated in 
e+e~ collision. Multiparticle exclusive amplitude Se^. .em{ki, . km), describing the 
production of m gluons : 

Se,..,Jk,,..,km) = (-ire-^nMp,(fc,)e-^; (7) 



i=l 



is a product of the m emission factors : 



MpXk.) = gs^^-^^QpAk^Gp-, (8) 
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where: 

gs = A/47ra;, 

n is the number of gluons emitted of quark (antiquark), 

ki — {oJii kj) denotes the 4-momentum of the ith soft gluon, 

Cj = ep\ J = 0, . . . , 3 describes its polarization, 

Pi is the 4-momentum of the parent of the ith gluon, 

Gp^ denotes the color factor for the given vertex of the tree diagram D, 



and can be conveniently represented as a tree diagram (see Fig. 1). Gluon emissions 
are not independent. Radiated particles have the memory of their parton parent, 
and of the previous parton splitting off its parent line. This dependence restricts the 
phase space of the produced gluon, and it is included in the form of a generalized 
step function Qp.{ki) : 

ep(k) := = a; < P°, ^kP < 0, ojO^p > Qo}; (9) 

where P is the momentum of the parent of a given parton k, 6 denotes the emission 
angle of the previous parton splitting on the P— line, and Qo is a cut-off parameter. 
Virtual corrections appear as a radiation Sudakov factor e , where : 

W{P) = J £k {A*p{k)Ap{k))^e) (10) 

denotes the total probability of emission of a gluon from the parent P, averaged over 
physical transverse polarizations e^,e^. Ap{k) is given by : 

It should also be emphasized that produced gluons are real (on-mass-shell) particles, 
so the energy cuk of a gluon of momentum k can be approximated as : 

a;fc =1 k I . (12) 

Summing of the color factors G over the color indices gives the result : 

/-Y /-Y* _ J Cp dla Pi — P , s 

^pAp. - I Cv dla Pi ^ P; ^^^^ 

where P denotes the 4-momentum of the quark (antiquark) which initializes the 
gluonic cascade. 
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In DLA different tree diagrams come from different non-overlapping kinematic re- 
gions, and do not interfere. Therefore, to calculate exclusive and inclusive multigluon 
densities it is enough to sum up incoherently the squares of amplitudes (0). Hence 
one obtains the exclusive density pff (ki, . ., k^) in the form : 

m 

p^^(ki,..,k„) =e-(^) En(^K(^^)^P.(^^))(eOe~'"^'^^ (14) 

D i=l 

parametrized by the momentum P of the quark (antiquark) which initializes the 
cascade. Multigluon inclusive density pp (ki, . .,km) follows from as : 

p5?(ki, . ., k^) = E / [dkU+,...n p^^Xki, . ., k^. (15) 

{n-m)\ J 

Introducing the method of the generating functional (GF) (see and references 
therein) allows us to perform the summation over diagrams in ([T^ and (|l^) in 



a very convenient way. While constructing the generating functional Zp[u\, one 
applies explicitly the self similarity property of the gluonic cascade. As a final result 
one obtains the recursive master equation in the form : 

Zp[u] = e-"'^-^) exp{ I d^k {A*p{k)Ap{k))(^,) u{k)Zk[u]). (16) 



It can be proved |^ that equation (|1^) reproduces contributions of all tree diagrams 
D, and allows us to express multigluon densities pp (ki, . .,km) and pp (ki, . .,km) 
as : 

p^^'(ki, . ., k^) = ^— i^Zp ||„=o}, (17) 

in 

Pp (ki, .., km) = — Zp\iu=n] (18) 

OUi . . . OUm 

where Ui denotes the probing function u{ki) and the functional derivative denotes 
sJkT) respectively. 

We would like to emphasize the simplicity of the GF approach. Especially for 
inclusive densities the method allows us to skip the complicated summation proce- 
dure, and to express the required distribution in a simple, compact form (for details 
see PI). 



2.3 Density matrix in the DLA formalism 

The DLA formalism in momentum space gives a good description of the structure of 
the gluonic cascade. The GF scheme suggests clearly how to construct multigluon 
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densities, and allows us to investigate their properties in a simple way. Since the 
method works so well for multiparticle distributions, one can expect to apply it 
successfully for other multiparticle observables. Below we briefly discuss the calcu- 
lation of the multiparticle density matrix in the framework of double logarithmic 
approximation, recalling our main results from Ref. 0. 

First we derive the general expression for the exclusive and inclusive density 
matrices dp (k'^, . ., k'^; ki, . ., k^) and (ip (k'^, . ., kj^; ki, . ., k^)- The task looks quite 
complicated because in this case the interference between different diagrams in (Q), 
generally does not vanish. Let us define two functionals, Zp[u] and Zp,[w], which 
generate the sum of all tree amplitudes and the sum of their complex conjugates 
respectively : 

Zp[u] = e-'"^^^/^exp{J d^kAp{k)u{k)Zk[u]), 

Z*p,[s] = e-'^^^'^/^expiJ d^kA*p,{k)s{k)Zl[s]). (19) 

The multigluon density matrices can be then expressed as (see Appendix A) : 
dp (k]^, . ., k^; ki, . ., k^) = 

= r T r Zp[u]Z*p,[s] |{^,=^ = 0},P=P' , (20) 

dSi> . . .dSm' OUi . . .dUm 



= Ssr...6s^,5u,...5uJ'^''^^*^'^'^ I{-^.^o|,p^p' • (21) 

Eq. ( pip is, in fact, a complicated integral equation. The difficulty of the diagram 
summation does not disappear there. It is only hidden in the compact form of (pT]). 
However, we remember that in the above formula we have taken into account inter- 
ference between all different graphs D and D'. And detailed analysis gives the result 
that in DLA not all the diagrams mix up: one can distinguish some interference 
classes. However, at the general level we did not succeed in formulation of such a 
GF which would take this fact into account. 

Nevertheless, we do not need the most general form of the density matrix. We 
are interested in its behaviour when the differences of momenta 
I ki — k[ I , • • • , I km — k'^ I are small, since we expect that large momentum 
differences will not contribute to Fourier transforms W%- It can be shown that 
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within this hmit interferences between different diagrams vanish, and one sums up 
only "squared" contributions from identical graphs. In fact, from the analysis of 
the diagrams contributing to single particle density matrix d'^p{k[,ki), it has been 
possible to prove (see Appendix B) that interference between different diagrams 
appears only if either energies or emission angles of produced particles are strongly 
ordered : Ui ^ Ui' {ui -C uJi') or 9ip ^ Oiip {6ip ^ Oyp). This statement can 
be generalized for any m-particle density matrix. If we have m particles, and (m-1) 
ones among them are "close" to each other, e.g. ki = k[, . . . , k^-i = ^m-D then 
the interference of the different diagrams will take place only if either the energies 
or the angles of km and k'^ are strongly ordered. 

Hence in our approximation we exclude mixing up different diagrams, and sum 
up only the squared contributions from identical ones. The exclusive and inclusive 
density matrices then take the simpler form : 

d-p (k]^, . ., k^; ki, . ., k^) = 

m 

= E {Se,..eJk[, . ., k'JS,,,,,Jk„ . ., U)(ei..e„), (22) 

D i=l 

dp (k]^, . ., k^; ki, . ., km) = 

cx) -| „ m n 

n=m \"' '"')■ D •' i=l j=m+l 

X ('S'ei. .e„(^l5 • • • ; ^m' • • • ; kn)Se^, .e„(^l; • ^n))(ei. .e„)- (23) 

The summation in (|23|), (p3D over D can be easily performed using the generating 
functional which reproduces only contributions of identical tree diagrams. In I 
proposed a master equation for such a generating functional (GF) Zpip[u{k', k)]. It 
accounts only for the contributions of "diagonal" diagrams (for proof see Appendix 
C), and generates the exclusive and inclusive density matrix as : 

c^p (k'l, . ., k^; ki, . ., km) = t t ^p'p |{«=o},p=p', (24) 

c^p (ki, . •, k^; ki, . ., km) = t 7. Zpip | ju^^an/.n} p=p/ . (25) 

The explicit form of its master equation, and the other formulae needed for further 
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calculations are presented below. The GF reads : 

OO 1 

Zp,p[n] = e-^(^' ^) ^ - / [dk'],,„n[dk],,„nu{k[,h)..u{k'^,K) 

x(A^,(A;;)Ap(A;0)(e,)--(^P'(^:)^p(^n))(e„) 

X^fc'ifclM- ■Zk'^kr.lu] Pv,..,n';l,..,n] (26) 

where the function Pi',. provides the requested parallel angular ordering ( see 
Fig. 5) for n particles in the form : 



Pl',..,n';l,..,n = ^ Q{9y_^pi > . . > dk'.^P') 'd{dki^P > ■ ■ > dki„p)] (27) 

(il,. .,in)(iP£Tm(l,. .,n) 

the product of the single particle amplitudes ( A*p,{k')Ap{k))(^e) averaged over gluon 
polarizations (see Appendix D) reads : 

{A*p,{k')Ap{k))(^e) = Ap>pik',k) = 

= If G*p,Gp \ I Qpik')Qp{ky, (28) 

ZTT VAuj'-^UJ'^ tipkup'k' 

and the radiation factor is given by : 

WiP',P) = ^in±^, (29) 

The other notations are the same as in section 2. 2. 

Generating functional ( ^6]) describes the sequence of "emissions" : 
{A*p,{k')Ap{k)) (^e) of two particles k, k' from two parents P, P'. If the profile function 
u{l',l) is equal to and 6^{l' — I) respectively, then the functional Zp'p takes the 
form : 

Zp,p[u = 0]= e-^(^'^), 
Zpp[u = 6\l' - 1)] = 1. (30) 

For the profile function of the form u{l', I) = v{l) 6^{l' — /) and P = P' our Zpip[u\ 



reduces to functional (0), as expected : 

Zpp{u{l', I) = v{l) 6\l' - I)] = Zp[v{l)]. (31) 
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2.4 Single particle density matrix. Leading terms 

From relations (pSj), ( [26| ) one obtains a simple integral equation for the inclusive 
single particle density matrix : 

d'^{k'- k) = J dhApp{s, s)d'^{k'- k) + 

+ fp{k\k)App{k\k)Zy^k[u = 5\l' -% (32) 

where fp{k', k) equals 

fp{k',k) = e-^'sC.\^-''^-'-''-^y (33) 
Introducing the notation : 

gp{k\ k) = fp{k\ k) App{k\ k) Zy^k\u = 5\l' - /)]; (34) 
we may write the symbolical solution d^p{k', k) of (R^ in the form : 



oo „ 

d'^{k',k) = Y. d^si. .d^SnApp{si,si). .As^_^s^_^{sn,Sn) gs,Ak',k). (35) 

n=0 •' 



We present Eq. ( pq) as a final result obtained in Ref. 0. The exact solution 
of d^p{k',k) and its detailed properties will be discussed below. We may replace 
App[s, s) in (RH) with its explicit form taken from Eq. 



„=o J Si J ZnOp^^ J Sn J 27r6^„_i„ 



[Sr. 



xgs^k'.k); (36) 

where / dflsi denotes the angular integration over the direction of vector Si, Bjj = 
Qsi,Sj and b = g'sCp- Factor gs{k' , k) defined by expression (P^), reads : 

X e "30 Qo ' e Qo Qo ' 

oo . ^ 111 

X e(0fci > . . > 0fc„)e(^fc.i > . . > Ok'n) 

X 6^,(1).. 6^,(^)6^(1).. efc(n); (37) 
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where h' = QgCv- Rewriting its phase space restrictions in terms of Bessel functions 
(see Appendix E) one arrives at : 



X 



_6|ln2%£-ln2%i££| _6'/2{ln2^^+ln2^^) 
' Qo Qo ' p Qo Qo ' 



oo n / finin{k,k') 



/ roo /•mm 
n=0 i=l ^-^0 Jmax{-g 



^( ,. . . .. ) 



fc,i-l{fc',i-l) ' 



fSk,i-l fSk',i-l / \N 1 

X / ^ dekiMxiOki) / dek'iJo{xi9k'i)) ■ (3^ 



The term which dominates in ( p7|) in the double logarithmic perturbative ex- 
pansion : 6 — > 0, P — > oo, 61n^ ^ = const, has the explicit form (for details see 
Appendix F ) : 



One may check (see Appendix F ) that this term iterated in (|36D produces the DLA 
leading contribution to the density matrix d^p{k',k). Namely, if we introduce the 
notation : 

.?'"(*'.^ . |(2.r/^/||-e,(.o../ff/5j^e_(.„) 

X9£>(fc',i); (40) 
then there is a relation : 

(r^{k\k) < d'p'^^\k',ky, (41) 

and for k = k' in quasi- diagonal limit : 

(fp{k',k) = d'p^^\k',k). (42) 

Both ( ^11) and ( ^2|) hold in the perturbative limit, as well. Hence we may write 
finally that : 

n=0 J Si J ZTTtlp-^ J Sn J Zntln-l,n 

xgW{k',k). (43) 
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The above result simplifies significantly the single particle density matrix approach 
0. We apply it in section 4. 2. so as to improve calculation of the single particle 
density in mixed coordinates. 

3 DLA in configuration space 

3.1 Fourier transform in transverse coordinates 

The DLA soft parton cascade starts from the initial parton of momentum P. Mo- 
menta of all particles produced from parton P refer to the P direction: they depend 
on the transverse and longitudinal momenta kx and ki (see e. g. [T^) taken 
with respect to the P axis. 

This dependence influences the structure of the cascade in configuration space. 
For the single particle amplitude Eq. (^) takes the form : 

Sp{x) = J d^ke'^^'-'^^'Spik) = I dkLd^kTC^'^'-'^'^Spik); (44) 

where S'p(k) (S'p(x) ) denotes the amplitude to produce a single particle with the 
momentum k ( at the space-time coordinate x ) from the initial particle of momen- 
tum P, and the 3-dimensional product kx reads : 

kx = /clXl + kx xt; (45) 

where indices L and T denote longitudinal and transverse components of momenta 
and coordinates with respect to the P axis. 

Since DLA amplitude Sp = Sp{u!,9pk) (0) depends explicitly on u =\ k \= k 
(0) and 9pk, one should express it as a function of kj^ and k-p. To proceed, let us 
recall that the simple form of Sp = Sp{k,6pk) has been obtained from the original 
QCD expression using the approximation of small emission angles 6pk 1 0. In 
this approximation one retains only the leading term in 6pk. Hence, restricting 
ourselves to the leading terms in 6pk, we interprete product kOpk as the transverse 
momentum kr and k as the longitudinal momentum ki. The correctness of the 
substitution k ^ ki for e**^* in (|4^) is not so evident, since it should be done in the 
argument of Fourier transform. However, if we restrict ourselves to finite, small time 
intervals, i. e. t -C (P^^)~^, it also works with a high level of accuracy. To see this, 
let us analyse the dependence of /c on /c^ and fc-r- Following (0), it can be expressed 
as : 

{uj=)k = \lkl + kl. (46) 

For intrajet cascades there is a relation kx <^ ki. Restricting ourselves to the 
leading term k ^ ki, we ignore the contributions of the order of t^. In DLA the 



12 



ratio has a limiting maximum value (P' denoting the emission angle of the previous 

parton splitting off the P-line (compare (^)). Maximal energy carried by the gluon 
equals P. Since that the approximation k ^ implies, that one neglects terms of 
the order of P6^. 

Using the leading term approximation, Eq. (|4^) can be rewritten as the product 
of 2 separate Fourier transforms: 1-dimensional FT of [ki, t — xi), and 2-dimensional 
FT of (kx, — Xt) in the form : 

Sp{xL,^T,t) = J dkLt'^^'-'-''^^ J £kTe-'^^''^Sp{k ^ kLM)- (47) 

Eq. ( ^71 ) gives us a good tool to investigate the space-time structure of DLA. We 
may now concentrate on particle distribution in transverse coordinates, which is 
physically the most interesting case. For the sake of simplicity, let us consider 
for t = the single particle amplitude 5'p(fc,XT) in mixed transverse space and 
longitudinal momentum coordinates [|13[ defined as : 



Sp{kL ^ k, XT, 0) = y d^kr e~'^^''^ Sp{kL ^ k, kx). (48) 

For the single particle exclusive distribution pp(/c, xt,0), defined as the amplitude 
square : 

p^P^(fc,XT,0) =|5p(fc,XT,0) |2; (49) 
one obtains simple expression in the form : 

p^^'(fc,XT,0) = j^Jd^kTd^k'^e-'^''^'''''r'>^^Sp{k,kT)SUk,K) 

= [ d^kTd^k'e-'^^^-^''^^''^d'j^(k,K;k,kT). (50) 

[271)^ J 

Relation ( pH) holds also for inclusive single particle density (for proof see Appendix 
G): 

p*p"(fc,XT,0) = ^y'd2M'fc^e-^(''--'^T)xTrf-(fc^k'^.fc^k^). (51) 

and can be easily generalized for any multiparticle distributions. In this paper we 
study the inclusive single particle density in mixed coordinates pp (fc, xt, 0) for the 
QCD gluonic cascade in DLA approximation. 
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4 Inclusive single particle density in configuration 



space 
4.1 Physics 

Momentum and configuration space description of a particle source are related to 
eacli otlier by tlie Fourier transform. Tfierefore one expects to extract from the 
observables in momentum space some qualitative information about their behaviour 
in configuration space. This observation applies, of course, to single particle density. 
In the DLA approach the inclusive single particle density Pp (k) [jlT| reads : 

= Im'' (52) 

where h = QsCf- For constant k the Pp (k) is concentrated around P direction (see 
Fig. 6)( skipping for now the Qp{k) limitations ). We expect, that the single particle 
distribution p^p{t,x.) in configuration space : 

P'p"(A;,xt,0) = J d^kTd^k'Te-'^''^-''''r^''^d^{k,k'T.;k,kTy, (53) 

will be concentrated around the P direction, as well. However, it has to obey the 
uncertainty principle. Since there are cut-offs of form (^ in density matrix, its 
Fourier transform will contain some oscillations, due to the restricted integration 
region. Moreover, relation (^ implies, that both densities Pp(k) and pp(A;, xx,0) 
have the same normalization, if integrated over d'^k and dk d^XT respectively. 

Remembering the above remarks, we propose now a technique for deriving the 
explicit form of inclusive single particle density in configuration space. 



4.2 Single particle density 

We may calculate the single particle density pp(fc,XT,0) as a Fourier transform of 
DLA density matrix (^). Let us make the transverse Fourier transform of both 
sides of Eq. (||). We obtain : 

P^"(A;,Xt(p),0) = -i^ / rf2A;T(P)t^'A;^(p)e-*('^-(^'-'^np))'^np)4-(fc,k^r;A;,kT) 
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x£H^,k'T(s);A;,kT(s)); (54) 

where indices T(P),T(S) correspond to the reference frame with the 2;— axis placed 
along the P {S) direction. The term of the lowest order {n = 0) takes then the 
form : 

Pp ^^^(A;,xt(p),0) = 
2h O 

" 2^^®^T < ^ < j^ddpMxTk9pk)f. (55) 



The numerical plot of ( p5|) is presented in Fig. 7. There is a limiting maximum value 
for xt = 0, as expected, and for large xt the function has a power decrease with 
the best fit exponent x^^"^. The result confirms our intuitive analysis from section 
4. 1. The term Pp *^^'*(A;, xt(p), 0) is concentrated around the P direction, however 

stronger than p^p^^\k, kx) from Eq. (|5^) . 

Derivation of terms of an arbitrary order in expansion (^) is more complicated. 
To proceed, let us analyse the last part of (^). It reads : 

/ — / o 1?" Q«-i (gn) 9n^ (k, k^(s) ; fc, kT(s) ) = 



= VirJ—j W^§^ e„_,(..)-i— e„(A:')e„(A:); (56) 

where the index (i) refers to the vector Sj. Taking into account phase space restric- 
tion and dominating contribution of ^-^^ — (see Appendix E) expression (E^) can be 
approximated as : 
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= ^l^(^)^ -R 0n-l(fc')0n-l(A:) X 

X / dXnXnJoiXnOkk') L dOnkM^nOnk) d9 nk' J q{x nO nk') ■ (57) 

Scheme (^) can be repeated iteratively in (^31). After some transformations one 
obtains finally (for details see Appendix E, H ) : 

p'p"(fc,XT(P),0)= (58) 

= 7^^®(~^ <k < P) L dOpk I dOpk' — / dLpk'kMkOk'kXx) 
Ztt a Ztc Jo 



P ( r°° r^Pk rffpk' 1 

X Jo(^ 86 ln(— ) |y dxxJo{x9f:k') j^^^ daJo{xa) j^^ da'jQ{xa')y 



— kQ[— <k<P) dOpk L dOpk' — / -rrz ^^MkOk'kXTj 

2vr 9 J9sL J9si 2-k Jo A{epk,9pk',9k'k) 



P fSpk ftlpk' 

x/o(a 86 ln( — ) { / da da' — - 

^ yoo 2nA{a,a',9, 



k'k) 



where 9k'k denotes the angle between momenta k', k, which in the reference frame 
with the z-axis placed along P direction takes the simple form : 



9k'k = V 9lk' + 0pk-'^ dpk'dpk cos{ipk'k)] (59) 

(fik'k denotes the azimuthal angle between kx and ki^, and A(a, 6, c) denotes the area 
of the triangle with sides a, 6, c. The numerical plot of (pB| ) is presented in Fig. 8 . 
There is a limiting maximum value for = 0, as expected, and for large xt the 
function has a power decrease with the best fit exponent x^'^'^^. The result confirms 
our intuitive analysis from section 4. 1. The term pp^^"* (A;, Xt(p), 0) is concentrated 
around the P direction, stronger than p^p^^\k,\<.T^) from Eq. (|5^) . The function 
oscillates around its power-law profile. As already mentioned in section 4. 1. , this 
effect is due to the sharp cut-offs of form (|^) which restrict the integration space. 

4.3 Properties 

Let us analyse the properties of ( pSj ) in detail. Comparing ( |55D with the single 
particle density in momentum space (|5^), we get the same dependence on the energy 
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k, as expected. The complicated emission term (see Fig. 9) : 

26 In 1^— j |y dxxJo{x6kk') J^^^ daJo{xa) j^^ da' jQ{xa') j = 

= [ d's'^Assik,k'T;k,kT) (60) 

corresponds to the logarithm In ^ of (|5^) . In fact, after integration of Pp'{k, Xt{p) , 0) 
over d'^XT one obtains the same result as for (^) integrated over d'^kj' (see Appendix 
I). The density p^p{k, xt(p), 0) is positively defined, as well. Using the Bessel function 
identities from Appendix H, Eq. (E^ can be rewritten in the form : 



Pp(A;,xt(p),0)= (61) 

= /fce(^<fc<P) E7^1n"(^)(/ dx,x,.. dx„x„) Y: 

e ^oiniy k Jo Jo m,..^=-oo 

( f0pk ]^ 

X < L dOpk Jrm{Xi9pk). -JmA^nOpk) Jmi + ...+m„{xTk9pk) [[ L daiJo^Xiai) ^ ; 



which implies positive definitness. For xt = Pp( A;, xt(p), 0) reaches its maximal 
value, as expected. Diagramatically, formula (^) represents the chain of indepen- 
dent emission factors (|60|) transformed onto the transverse xt plane by the Bessel 
factor of primary emission from the parent P, namely Jo(| — k'y | Xt) (see Fig. 
10). For k = P expression (^) reduces to (^) with the dominance of the primary 
emission, as expected. 



5 Summary 

We considered the QCD parton cascade created in e~^e~ collision in the double loga- 
rithmic approximation. Using the density matrix (DM) formalism 0, and restricting 
ourselves to the terms leading in the double logarithmic (DL) perturbative expan- 
sion for the quasi- diagonal limit k' = k, we derived the explicit form of the inclusive 
single particle density matrix d^p{k\k) and single particle density pp (fc, xt(p), 0) 
(see Fig. 8) in mixed coordinates (Jzl ~ k,ii.T)- The gluon density pp (fc, xt(p), 0) 
fulfills important physical requirements such as positive definitness and proper nor- 
malization. It is concentrated around the P direction, and shows the power law 
profile for large x^. Due to the cut-offs of the type @j which restrict kinematic 
regions, the density oscillates around its power law profile. 
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The above results give a positive outlook for the future. The simplified technique 
for calculating multiparticle observables in DLA for the constant 0:5 is ready. It may 
allow one to investigate the structure of the QCD cascade in a very comprehensive 
way: the exact form of a Wigner function already would give us clear experimental 
predictions. 

In order to improve the approach one needs to include the momentum dependence 
of as- The exact analysis of applicability of the quasi- diagonal approximation 
k' = k would be required, as well. There is also a problem as to what extent the 
kinematical constraints characteristic for a given QCD approach ( in DLA they are 
of the form ) influence the results obtained in the approach. In other words, 
there is a problem how to separate the effects of the kinematical restriction from 
the dynamical content of QCD space. The explanation would require performing an 
analysis similar to the presented above for LLA and MLLA schemes, as well. The 
MLLA approximation would be of special interest since it would incorporate the 
energy conservation into the cascade. 



6 Appendix A 

The functional Zp[u\ and 2'p/[s] defined in (|19]) generate the sum of amphtudes 
and their complex conjugates over all possible tree diagrams respectively. It can be 
seen when one rewrites the master equation for Zp in the form of the diagram series 



(see Fig. 2). From that construction follows exclusive density matrix (20) : 



dp (k'l, . ., k^; ki, . ., k„) = Z*p[s\ |{^=o} t 7 — Zp[u] |{„=o} • (62) 

OSii...OSjn' OUi...OUm 

Then the inclusive density matrix calculated from (|^) looks like this : 

001/ A™ 

d-?(k;....k-k,....k.„) = i:- (/A--] ;^_^zM.| |,..„, 

X 7 7 — Zp[u] \{u=o} ; (63) 

OUi . . . OUm 



and from the identity for the product of any two functionals F, F' : 



(64) 



follows {M 
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7 Appendix B 



In DLA there are four different tree grapfis Ma, Mh, Mc, M^, describing the produc- 
tion of two gluons. They are defined on non-overlapping kinematic regions 
(see Fig. 3 ). Emitted gluons are either angular (AO) or energy ordered (EO). 

Let us consider all the diagrams contributing to the single particle density matrix 
d'p{k[,ki). From the (AO) and (EO) it follows that the interference between any 
two different graphs will appear only if either energies ui , uj[ or emission angles 
Oip, 9i'p of produced gluons are strongly ordered (Fig. 4). 

This statement can be generalized for any m-particle density matrix by induction. 
If we have m particles, and (m-1) ones among them are "close" to each other, i. e. 
ki = k[, . . . , km-i — ^m-i; then the interference of the different diagrams will take 
place only if either energies or angles of k^ and k'^ are strongly ordered. 

8 Appendix C 

Z p/ p[u{k' , k)] defined in (^) generates correct exclusive density matrices (0). The 
proof of this statement follows from the representation of the master equation for 
Zp'p as a diagram series (see Fig. 5). The series reproduces all squared contributions 
of identical tree diagrams, and excludes interference of the different ones (function 

Pl',...,n';l,...,n)- 

The explicit form of inclusive density (^) follows from formula (0). Substituting 
(P^) into (^) one obtains : 

(^15 • -5 ^m) ki, . ., km) = 

= E ^ / d'kd'k'6\k' - k)- Zp,p[u] |{„=o}; 

(65) 

which represents the functional Zpip[u\ expanded around "null" {u = 0} in the 
"point" u = 5^ir-l). 



9 Appendix D 

We want to average relation (pSD over physical polarizations. The exact expression 
to be summed over polarizations looks like this : 

(e-P)(e'.P'). 



{k-P){k' -P') 



(66) 
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where e, e' are polarizations of the same gluon identical in the limit k' = k. The 
gauge fixing we apply in the approach allows us to neglect contributions to ( |66D 
coming from the " nonphysical" polarizations e° and e^. Furthermore, it requests 
the time components of the physical polarizations e^,e^ to be equal to 0. Space 
components of e^, can be then constructed in the form : 

P X k 



= (67) 



and for e' respectively 



e'2 



1 p 


X 


k 1 ' 




X 


k 


1 e2 


X 


kl' 


P' 


X 


k' 


1 P' 


X 


k' 1 




X 


k' 


1 e'2 


X 


k' 1 



= (68) 



Summing over these 2 polarizations and expanding scalar products in (pq) , one 
obtains finally : 

(e' ■ P)(e" ■ P') ^ J_ ^_ 
(k-P)(k'-P') uiuj' 9pt9p.»/ ^ ' 

The result can be easily confirmed for any two physical polarizations e^, placed 
on the plane e^e^ : 

= cos(^ + smy9 , 

= — sirup + cos(f e^; (70) 
and for e'^, e'^ placed on the plane e'^e'^ with the phases ip' = (p respectively. 



10 Appendix E 

To transform (|56|) into (|57|) we will apply the pole approximation method ([|r^ and 
references therein). The phase space restrictions 9„_i(?2)9„(A;)6„(A;') in (^), in 
particular the angular ordering (AO) 9n~i,n ^ On,k, 6n~i,n ^ Qn,k' make the the 
term ^ ^ be dominating in (p^. Applying the Bessel function identity : 

1 r dOnk dOnk' ^nk ^nk' 



2tX 2iT J A{9nk,dnk',dkk') 

dxxJo{x6k'k) / dOnkOnk Jo{x6nk) I dOnk'^nk' Jo{x6nk')', (71) 
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where Oyk denotes the relative angle between vectors k, k', and A{6nk,(^nk',(^kk') 
equals the area of triangle with sides 9k'k, Gk'm Gkm we may rewrite expression (|56D 
in the form : 



271 J S„. J 27l9l,^ " "^'"^ Onk'Onk 



2b r dsr, r°° 



27ik^ 
1 



Joi^nOkk') / dOnkJoi^nOnk) / d9nk' Joi^nOnk') 
J Sn Jo J J 

■e„_iHe„(A;')e„(fc). (72) 

Because of angular ordering the angle 6n-i,n practically does not change while in- 
tegrating over angles Onk,Onk', and can be successfully approximated by the angle 
&n-i,k (^n-i,fc')- Applying all the integration restrictions explicitly, we arrive finally 
at the expression : 



2Tck^ J Sn J 2'n9i_i Onk'^nk 



^ e„_i(fc')0n-i(A;) r^^x (73) 



2Tlk^ 9n-l,k'0n-l,k Jk S 

X dXnXnJo\XnOkk') J d6nkJo{Xn6nk) J d6nk'Jo{XnO' 

which after trivial integration over s„ gives identity (|^. 



nk' 



11 Appendix F 

We shall find such terms in expansion (|36|) which dominate in the double logarithmic 
perturbative limit of d^p{k',k), i. e. when b ^ (6 oc as) and P oo, so that 
61n ^ remains constant, and generate double logarithmic corrections to the cross 
section. First let us introduce the following notation : 

oo 

d^^{k',k) = (74) 

n=0 

where coefficients a„ describe the nth order iteration of gs{k', k) (|36D in the form : 
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The term of the series ( |75D is equal to gp{k', k). Expanding it in the powers of 
couphng constant b, one obtains the term of the first (lowest) order of b in the form : 

Denoting the other terms by the symbol Oq^^ = g^p^\k', k), the term ao from the 
series (|75D can be rewritten as : 



ao = + ai^'l (77) 

We have checked that Og^'' dominates in (^). To see this, one should rewrite ao as : 

ao = gp{k\k) = aj,') + af.^'^ (78) 
= a^^f{k\k-P)- (79) 

where the exact form of /(A;', k\ P) follows from (0). For 6 > and P finite the 
f{k',k;P) vs I k — k' I is a gauss-like function with the maximum equal to 1 for 

PT 

k = k' and a non-zero minimal value. In the DL perturbative limit : f{k', k; P) 1. 
Consequently, for any k, k', P : 

ao < al^\ (80) 



a. 



(1) 





1. (81) 



Since ao^"* dominates in ao, we expect that the iterations of ao^"* will generate the 
leading contributions to the dp (fc', k). Let us introduce the notation : 



oo 



d'-'''\k',k) = Y.<\ (82) 

n=0 

oo 

d'^^^>'\k\k) = (83) 



n=0 



where a^-*^-* and a^'^ represent the result of nth iteration (^) of ao^^"* and ao^'* re- 
spectively. From (pOD immediately follows the relation : 



an = a« + a(fi) < a«. (84) 
Since a„ > 0, for the density matrix one obtains : 

< (rp{k',k) = d'^'^^\k\k) + df^''^\k',k) < d';i'^'\k',k). (85) 
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In the limit | k — k' | ^ the density matrix d^p(k', k) — > (fp'^^\k', k), and there- 
fore : 

dfjk', k) PT ^ ^gg^ 
dp (^'7 ^) 

For other ( larger ) values of | k — k' | the contribution of (fp'^^\k',k) generally 
need not to dominate. However, taking into account the fact that the density ma- 
trix approach works only for the quasi- diagonal limit, we may identify the quasi- 
diagonal region with the region of the dominance of (fp'^^\k', k). Furthermore, since 
(fp'^''^\k,k) = 0, the (fp'^^\k',k) generates all DL corrections to the cross section. 
Hence we finally arrive at the relation (^). 

12 Appendix G 

The particle density pp(/c,xr,0) can be represented as : 

Pp(fc,Xr,0) = J2 ( _ i\l / dk2(fx2T---dknd^XnTP''pik,XT,k2,X2T,kn,XnT)- 
n=l 

(87) 

Representing the exclusive particle density as the square of multiparticle amplitudes, 
expressing the configuration space amplitudes as Fourier transforms of amplitudes in 
momentum space, and performing integration over d'^X2T ■ ■ ■ d^XnT finally one arrives 
at Eq. (^. 

13 Appendix H 

For our purposes let us quote the following identities for Bessel functions (|T^ and 
references therein) : 

00 

Uxo^k) = E e''^^''~^'^J^{xOkp)Jm{xek'p), (88) 

m=— 00 
00 

gifcxcos^ _ ^ i'"e''"^J„(A;x), (89) 



00 







/ \ / i\ ^{'^ — a') 
dx xJo{x a)Jo{x a ) = ; (90) 



where 6k' k denotes the relative angle between momenta k' k, which in the reference 
frame with the z-axis placed along P direction takes form (^), and (fk'k = ~ f\ 
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where the angles y?, denote the azimuthal angles of vectors kx and k!^ on the 
transversal plane respectively. 

Now let us repeat the scheme Eq.(^)-^ Eq.(^) iteratively in (^3|). Then, after 
introducing the explicit form of Fourier transform we arrive at the expression : 

p^p"(fc,XT(P),0)= (91) 

If r /-Stt /'27r 



//•ZTT pZTV 
d{k9kp)k9kP difi ^^f^-ike^pXT cos ^+ikeypXT cos ^ 
Jo Jo 



(27r)2 

(2&)" l.n fP\ 26 1 



oo 



nl \kj 2Tck-^ 0pk9pk' 

X (^y^ dXiXiJoiXiOk'k)-- dXnXnJo{XnOk'k)j 



/ daiJo^xiai) / da[Jo{xiai>). . 1^ (ia„Jo(x„a„) / ' (ia„Jo(x„a„) ^ 

It looks still quite complicated, however one can notice that the term in brackets : 

/ daiJo{xiai) I rfa^ Jo(xiar). . / (ia„Jo(x„a„) / da^Jo{xnaJ > ; 

(92) 

which contains integration of 2n Bessel functions, simplifies significantly if one adds 
to ( |9TD some extra off-diagonal terms, which are equal to in the quasi-diagonal 
limit. Namely : 

1 r^Pk fSpy fSpk 



£)Ly^ 1 rfpk f"Pk' rfpk f^Pk' 

{} = — daiJo{xiai) L (ia^ Jo(xiar). . / (ia„Jo(x„a„) / c/a^ Jo(x„a^). 

(93) 

Then expression (^) takes the simpler form : 



p^"(fc,XT(P),0) = ^kQ{^<k<P)j^^x 

X /n dOkP C dOk'P r dip j'^^^'e-^^dupXTCos^+ikeypXTCos^' /g^^, 
J9a j9a Jq Jo 



X ^ In" 1^— j dxxJo{x9k'k) daJoixa) da'Jo{xa) 
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Let us perform the integration over the angles Replacing all the terms of the 

type Jo{x6k'k) by their expansions (0), expanding then Fourier transforms in terms 
of (|89D , and performing the integration over Lp, Lp' explicitly, we arrive at (|6ll) . 



However, the part of expression (^) containing Bessel functions can be rewritten 



as : 



n oo 



W X! JmX^i^Pk)Jmi{XiOpk') ] Jmi+. .+mn{^TkO pk) Jmi+. .+m„{xTkO pk') 
\i=l mi=—oo 



d(Pk'k\Y[M^i0k'k)jJoixTk6k'k); (95) 



1 /•2^ 
271 Jo 



where cpk'k denotes the relative angle between the vectors kj-, and 



Ok'k = y^pk' + ^Pk ~ 2 6'pfc/6'pfc cos(v3fc'fc), as usual. Substituting ( p5| ) into ( ]61| ) finally 
we obtain (158|). 



14 Appendix I 

Let us investigate normalization of single particle density ( pHD in real space, i. e. the 
integral / (i^XTPp(fc, x^-, 0). It reads : 

d^XTffp{k,XT,0) = 

2bkQ{— < k < P) N ^ dOpk j^^^ ^^"'^ dVk'kj dxTXTJoikOk^kXr) 



x/o(^ 86 ln(— ) |y dxxJQ{x6kk') daJo{xa) da' jQ{xa') j) . (96) 

The term JoikO^ikXT) is the only part of expression ( PB| ) depending on xt- Integral 
dxTXTJoikO^ikXT) can be performed, using identities (|88|) , (pO|) . It equals : 

, dxTXrJoikdk'kXT) = 27r 5(v9fc/fc)^^%^— (97) 
k-^OkP 



Substituting { p7\) into (|96D , and integrating over 6k' p and (pk'k, one arrives at : 

2&^,Qo , „^ /-^ dOpk 



J d'xTP^{k,^T,0) = ^©(^ < ^ < 



X 



Pfc 



x/o(^ 86 ln(— ) |y dxx j^^ daJo{xa) j^^ da'Jo{xa')j). (98) 
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Performing the integration over x in the argument of function Jg with the help of 
identity (pOD, finally one obtains : 

l^.r^(M ^ f < * < P)/; ^/o(/8Mn(f)M^)); (99) 

which equals the normalization of particle density in momentum space defined in 
(0): 

J d^T Pp (fc, XT, 0) = y d^kr Pp{k, kr). (100) 
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Figure captions 



Fig. 1 Feynman diagram for the production of m gluons in DLA, where P denotes 
the 4-momentum of the initial q{q) and ki denotes the 4-momentum of the ith pro- 
duced gluon. 



Fig. 2 Generating functional ([T9|) as a diagram series 



Fig. 3 Redefined kinematical regions in DLA. 

Fig. 4 Interference between different diagrams for d']¥{k[, ki). Remark: diagrams 
(a) and (b) are identical except of the position of ki {k[) leg. 

Fig. 5 Master equation for generating functional (p6D represented as a diagram 
series (for details see 0). Function -Pi',...,n';i,...,n introduces the parallel angular or- 
dering (AO). 



Fig. 6 Function pp(k) from Eq. (|52D vs kx =\ ky | in Qq units for parameters 
b = 0.25, P/Qo = 243, = 1, k/Qo = 128, chosen following W^. The power fit 



reads p'^ik) = (g;)"^'^^ * 0.00065 Qq - Plot (a) with the logarithmic scale for the 
vertical axis, plot (b) with the logarithmic scale for both vertical and horizontal axes. 

Fig. 7 Function p^p'^^\k,:Kt,0) from Eq. ( ^5]) vs. xt =\ | for parameters as in 



Fig. 6. The power fit reads pf{k, xt, 0) = {QoXt)'^-^^ * 4:.0E - 05 Qo- Plot (a) with 
the logarithmic scale for the vertical axis, plot (b) with the logarithmic scale for 
both vertical and horizontal axes. 

Fig. 8 Function pp(A;,Xf,0) from Eq. (p^) vs. Xt =| xt | for parameters as in 
Fig. 6. Plot (a) with the logarithmic scale for the vertical axis, plot (b) with the 
logarithmic scale for both vertical and horizontal axes. 
The power fit reads p'^(fc,xr,0) = (QoXt)'^-'^^ * 0.00018 Qo- 

Fig. 9 Diagramatic representation of the term 

/P\ ( f°° f^pk r^py 



26 In ^— j |y dxxJo{x6kk') J^^ daJo{xa) j^^ da'jQ^xa) 
dh — Assik,'k'T;k, kr); 
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The integration region lies in the overlap of conus Qs{k) and conus Qs{k')- 

Fig. 10 Diagramatic representation of p^p{k, xt, 0) from Eq. (^) as the chain of 
independent emission factors ( ]60| ) transformed onto the transverse xt plane by the 
Bessel factor of primary emission from the parent P, namely Jo(| — k'^ \ ^t)- 
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